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Abstract We investigate the existence and the global causal structure of plane 
symmetric spacetimes with weak regularity when the matter consists of an irro- 
tational perfect fluid with pressure equal to its mass-energy density. Our theory 
encompasses the class of W ' regular spacetimes whose metric coefficients have 
square-integrable first-order derivatives and whose curvature must be understood 
in the sense of distributions. We formulate the characteristic initial value problem 
with data posed on two null hypersurfaces intersecting along a two-plane. Relying 
on Newman-Penrose's formalism and expressing our weak regularity conditions 
in terms of the Newman-Penrose scalars, we arrive at a fully geometrical formu- 
lation in which, along each initial hypersurface, two scalar fields describing the 
incoming radiation must be prescribed in L 1 and W~ 12 , respectively. To analyze 
the future boundary of such a spacetime and identify its global causal structure, 
we introduce a gauge that reduces the Einstein equations to a coupled system of 
wave equations and ordinary differential equations for well-chosen unknowns. We 
prove that, within the weak regularity class under consideration and for generic 
initial data, a true spacetime singularity forms in finite proper time. Our formula- 
tion is robust enough so that propagating discontinuities in the curvature or in the 
matter variables do not prevent us from constructing a spacetime whose curvature 
generically blows-up on the future boundary. Earlier work on the problem studied 
here was restricted to sufficiently regular and vacuum spacetimes. 
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1 Introduction 

Our objective in this paper is to establish the existence of a large class of space- 
times with weak regularity, and to shed some light on the propagation and interac- 
tion of gravitational waves. Specifically, we consider plane symmetric spacetimes 
satisfying the Einstein field equations of general relativity when the matter con- 
sists of an irrotational perfect fluid, whose pressure equals its mass-energy density. 
This latter assumption implies that the sound speed coincides with the light speed 
and, therefore, discontinuities in the fluid variables may occur only along null 
hypersurfaces. 

Our aim is to construct a large class of plane symmetric spacetimes by for- 
mulating the characteristic initial value problem when data (with weak regularity) 
are prescribed on two null hypersurfaces intersecting along a two-plane. In addi- 
tion, we investigate the regularity and global causal structure of such spacetimes, 
and establish the strong version of Penrose's cosmic censorship conjecture for 
spacetimes with weak regularity. Earlier work on this problem was restricted to 
sufficiently regular and vacuum spacetimes. As far as the matter model under con- 
sideration is concerned, the initial value problem was tackled earlier, especially by 
Taub [ 39 1 (see also [ 38 1 ) who considered the initial value problem on a comoving 
hypersurface and derived an integral formula for sufficiently regular solutions. 

We are primarily concerned with the characteristic initial value problem and 
its well-posedness when, on each of the two hypersurfaces, two scalar fields rep- 
resenting the gravitational radiation and the mass-energy density of the fluid are 
prescribed. Our presentation, based in part on the Newman-Penrose (NP) formal- 
ism, yields a geometrical formulation for the characteristic initial value problem 
in which physically relevant (i.e. chart-invariant) quantities are clearly identified. 
With regard to earlier related works, we want to mention that the characteristic 
initial value problem, addressed in the present paper for matter spacetimes with 
symmetry and low regularity, was treated earlier within the class of sufficiently 
regular vacuum spacetimes by Friedrich ifTTIl . Stewart and Friedrich |f36l , Ren- 
dall |32), and Christodoulou lfl2l 

We now outline what is meant by "weak regularity" in this paper. First of all, 
we emphasize that throughout Sections 2.1 to 3.4 we assume the spacetime to be 
regular. In Section 3.1 we show that spacetimes of the symmetry class studied 
here are determined by the solution of singular wave equations (cf. (13.4b . below), 
for which the concept of weak solution is well-defined. We shall call spacetimes 
generated by these weak solutions "weakly regular". Equations of the form d3 .4b 
have been studied extensively and are know as Euler-Poisson-Darboux equations. 
We shall write down an integral representation of solutions of a certain charac- 
teristic initial value problem in terms of the Riemann-Green function, which is 
known explicitly, and the initial data appropriate to this problem, thus giving a 
representation of solutions both strong and weak solutions. 

Most importantly, our results assume weak regularity on the data and, there- 
fore, on the solutions, so that the Einstein equations must be expressed in the 
sense of distributions. Specifically, we introduce the class of W 1,2 regular space- 
times by requiring that the metric coefficients (in a certain gauge) belong to the 
Sobolev space W ' of functions with square-integrable first-order derivatives. The 
low regularity imposed on the metric still is sufficient for the Einstein equations to 
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make sense within the theory of distributions |23~1 . The Riemann curvature tensor 
then belongs to the dual space W~ 1 ' 2 consisting of distributions that are derivatives 
of square-integrable functions. 

More precisely, in view of the Einstein equations, the Ricci part of the curva- 
ture turns out to be solely in L 1 (that is, integrable) and may contain propagating 
jump discontinuities (along null hypersurfaces of the spacetime), while the Weyl 
part is solely a distribution in the negative Sobolev space W~ 1 ' 2 and, for instance, 
may contain Dirac masses supported on null hypersurfaces. For further results on 
the existence and qualitative properties of weakly regular spacetimes, we refer to 
LeFloch and Rendall lf24l . LeFloch and Smulevici l25l . and LeFloch and Stewart 
11261 (see also 0). 

The importance of spacetimes with very low regularity was first recognized by 
Christodoulou, who constructed (cf. [8|-[11|) spherically symmetric spacetimes 
with matter, provided a detailed description of their global structure, and estab- 
lished the strong version of the cosmic censorship conjecture (within the restricted 
class of spacetimes under consideration). In comparison, the plane symmetry as- 
sumption made in the present paper leads to a technically simpler set of partial 
differential equations but yet allows us to study gravitational waves (and their in- 
teractions) which do not arise in spherical symmetry. 

To investigate the global structure of the spacetimes under consideration, we 
introduce a gauge in which the reduced Einstein equations take the form of a sys- 
tem of wave equations with singular coefficients coupled to ordinary differential 
equations. The key equation is of the Euler-Poisson-Darboux type which we solve 
with the Riemann function technique. The novelty of our analysis lies in the weak 
regularity assumed by the solutions, as the initial data solely belong to the space 
W 1 ' 2 and solutions must be understood in the distribution sense. Uniqueness holds 
in the symmetry and weak regularity class under consideration. 

Within the low regularity class under consideration, we construct a develop- 
ment of a given initial data set and we study the possible formation of curvature 
singularities in this spacetime as one approaches its future boundary: we prove that 
generic characteristic data always lead to a curvature singularity. To this aim, we 
study the blow-up behavior of the Riemann function and derive the leading behav- 
ior of certain curvature scalars, which generically are found to, indeed, blow-up 
to infinity near the future boundary of the constructed spacetimes. For sufficiently 
regular and vacuum spacetimes, a related result was established earlier by Mon- 
crief [27 1 for Gowdy spacetimes with torus topology. Typical examples of this 
behavior (in the vacuum) were also constructed earlier explicitly by Khan and 
Penrose i2TTl and Aichelburg and Sexl Q. 

Interestingly, the spacetimes constructed in the present work can be interpreted 
as "colliding spacetimes" in which two plane gravitational waves propagating in 
a flat Minkowski background collide. Such spacetimes have been constructed ex- 
plicitly for special choices of initial data sets [2,6], while other authors inves- 
tigated general stability properties of colliding spacetimes [22 3 1 30, 37], even 
going beyond the plane-symmetric case II41II42I . Furthermore, the characteristic 
initial value problem was treated in [18,19] when the initial data set is sufficiently 
regular. More recently, the fomation of trapped surfaces was analyzed in H15II16I 
[40). All of these earlier works were concerned with regular initial data sets, and 
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our results in this paper generalize some of these results to the broad class of 
weakly regular spacetimes. 

An outline of the structure of this paper follows. In Section 2 we describe our 
matter model and write down the Einstein field equations for a regular spacetime. 
In the first part of Section 3 we set up a characteristic initial value problem where 
the evolution equations are two uncoupled linear second order hyperbolic equa- 
tions, and we obtain an explicit integral representation of their solutions in terms 
of appropriate initial data. There is a well-defined concept of weak solutions for 
these evolution equations and this leads in the second half of Section 3 to our def- 
inition of W 1,2 space-times, and Theoreml, the existence theorem. These results 
are local and depend on a particular choice of chart, and so in Section 4 we recast 
the results in chart-independent form using the Newman-Penrose (NP) notation. 
In Section 5 we obtain our global existence theorem and investigate the details 
of curvature blowup. In the final section we investigate jump conditions across a 
null hypersurface where the metric fails to be regular and study the propagation of 
curvature singularities. 



2 Einstein equations for plane symmetric spacetimes 

2.1 The Einstein field equations 

In this section we first present our assumptions and write down Einstein field 
equations for the geometrical set-up and matter model under consideration; we 
mostly follow the presentation in Tabensky and Taub 113 81 . Throughout Sections 
2.1 to 3.4, the metric and fluid variables are assumed to be regular solutions of 
the field equations we are deriving. Starting in Section 3.5 we consider weak solu- 
tions (defined in the sense of distributions). We are interested in plane symmetric 
spacetimes [j^,g) — or polarized Gowdy spacetimes — described by 



where the scalars a,b,c depend only upon the characteristic variables u = t — x, v = 
t +x. Throughout this paper we use the signature (+, — , — , — ), which is standard 
in the physics literature. Observe that, while the functions a and c are coordinate- 
dependent, the function b carries a geometric-meaning and e h represents the area 
element of the surfaces of symmetry described by the coordinates (y,z) and, if 
the coordinates on the orbits of symmetry are rescaled so that (y,z) *-> (ky,kz) for 
some k > 0, then e h h-> e b /k. 

The Einstein tensor G a p associated with the metric J2.lt , together with the 
energy-momentum tensor of the fluid T a p , must satisfy the field equations 



g = e 2a {dt 2 - dx 2 ) - e 2b {e 2c dy 2 + e- 2c dz 2 ) 
= e 2a dudv - e 2b (e 2c dy 2 + e" 2c <fe 2 ) , 



(2.1) 



G a p — K T a p, 



0,0=0,. ..,3, 



(2.2) 



where K > is a constant. It is straightforward (but tedious) to compute the 
Christoffel symbols and curvature coefficients associated with the metric J2.lt . 
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and to arrive at the following expression of the relevant components of the Ein- 
stein tensor: 

G 00 = 2 (2a u b u - b m - b\ - c 2 .) , 
Goi =2(b uv + 2b u b v ), 

G n =2(2a v b v -b vv -bl-c 2 v ), (2.3) 
G 22 = -4e- 2a+2b+2c (a uv + b uv + b u b v - b u c v - b v c u - c m + c u c v ) , 
G33 = -4e- 2fl+2fo - 2c (a uv + b m + b u b v + b u c v + b v c u + c uv + c u c v ) , 

where the subscripts on a, b and c represent partial derivatives with respect to u, v. 
Recall that indices are lowered or raised using the metric, for instance G a p : = 

gaygpBG 75 - 

We assume that the matter is irrotational null fluid, that is, an irrotational per- 
fect fluid whose pressure p is equal to its mass-energy density w, i.e. 

p = w. (2.4) 

It is described by the energy-momentum tensor (with our choice of signature) 

T a P = {w + p)u a uP -pg a P, (2.5) 

where u a denotes the 4-velocity vector of the fluid, normalized so that 

u a u a = 1. 

As we will see below, the condition p = w implies that the sound speed in the fluid 
coincides with the light speed, normalized to be 1 . 

The (contracted) Bianchi identities for the Einstein tensor, V a G a ^ = 0, in 
combination with the field equations (12.2b are equivalent to the Euler equations 

VaT"? =0. 

Under our symmetry assumptions and for the matter model under consideration, 
the Euler equations are equivalent to the following partial differential equations 

2wMU a uP +2wu a , a uP +2wu a ift- a -w M g afi =0, (2.6) 

where, for instance in u^- a , the subscript ;ce denotes the covariant derivative. 
On the one hand, we can multiply ( 12.6b by up and obtain the scalar equation 

2w M u a — 2wu a - a +2wu a Uj } u^- a — w, a u a = 0, 
which, in view of the identity upu^ — 0, simplifies into 

U a W^a+2wu a - a = 0. 
Assuming that the density is bounded away from zero and setting 

L = Jflogw, 
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we conclude that 

u a La+u a . a =0. (2.7) 

On the other hand, we can multiply ( 12.6b by the projection operator Hpy — 
gpy — upiiy and obtain the vector- valued equation 

or equivalently 

H a ^L a - u a uf, a = 0. (2.8) 

In addition, we assume the flow to be irrotational and we introduce a potential 
\f/ associated with the velocity, whose gradient is timelike \jf p y>P > 0: 

Vm (2.9) 



W.p W'P 



The fluid equations ( 12.7b and ( 12.81) then become 

= (2.10) 



w 1 / 2 1//'™ 



W.p V 4 



and 



(E-log^/Y^) p =k\l/ ; p, (2.11) 
respectively, where the scalar k is 



The second equation, ( 12.11b . states that the gradient of E — log ^Y' a ¥,a is 
parallel to the gradient of (//, so that the former can be expressed as F(\jf) for 
some function F. By replacing \ff by some function G(y/) if necessary we can 
always arrange that E — log ^\ff' a \ff, a = 0, in other words 

w = yr a w (2.12) 

This is the relativistic analogue of Bernoulli's law for irrotational flows in classi- 
cal fluid mechanics. It determines the mass-energy density algebraically, once we 
know the velocity of the fluid. 

Finally, the equation ( 12.10b determines the evolution of the remaining fluid 
variable, that is, the potential \jf. By using the short-hand notation \j/ a = Y,a, Wu — 
y/ u , etc., and in view of 

Wa = (y«,Vv,0,0), iff" =2e- 2fl (v^., V«,0,0), 

it follows that 

W a ¥a =4e" 2fl y«Vv, 

:(Vv,Vk>0,0), U a = =(Yu,Yv,0,0)- 



\/Wv 2^/vWf 
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Hence, the equation ( 12.121) becomes 

H> = 4e- 2Q VW„, (2.13) 
while ( 12.101) reduces to a wave equation^ for the potential: 

□ V = Vwv + b v y/„ + b u y/,, = 0. (2. 14) 

The latter is the essential matter equation to be investigated. 

Two main assumptions were used in our derivation: we needed that w remains 
positive, and that the solutions are sufficiently regular. We will show how to relax 
the first of these in Section 2.2. The regularity of solutions will be discussed in 
Section 3.5 and subsequent sections. 

We are now in a position to write down Einstein's field equations for the ge- 
ometry variables. The components of the tensor T a p = w (2u a up — g a p) are 

r 00 = 2v/„ 2 , 7bi = o, r n =2v/ 2 , 

T 22 = Ae- 2a+2b+2c ¥uVv , T 33 = 4e- 2a+2b ~ 2c ¥llVv . 

Returning to (12.2b and relying on the expressions ( 12.3b of the Einstein tensor, we 
arrive at the (evolution and constraint) equations for the metric coefficients a,b,c : 

Kyi, 

0, (2.15) 

Observe that the first three equations contain second-order derivatives of b, while 
the last two equations are equivalent to the system 

c uv + b u c v + b v c u = 0, (2. 16) 



a m - b u b v + c u c v = -K\ff u \if v , (2.17) 

which contain second-order derivatives of c and a, respectively. At this stage of the 
analysis, we observe that all of the equations under consideration are nonlinear and 
involve quadratic products of first-order derivatives of the metric coefficients. This 
completes the derivation of the Einstein equations in characteristic coordinates for 
plane symmetric spacetimes (^,g). 



2,ci u b lt b m[ b u c u — 
2a v b v -b vv -b 2 - c 2 = 
b uv +2b u b v = 

&uv b uv b u b v b u c v -\- b v c u 4- c uv c u c v — 
a m + b uv + b u b v + b u c v + b v c u + c m , + c u c v = 



This is a general fact for null fluids, irrespective of our symmetry assumption. 
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2.2 Physical meaning of the matter model 

The reader will have noticed that our derivation above assumed that w defined 
in (12.12b remains positive. This condition may be imposed initially on the given 
data but, in general, will fail after some finite time. When this happens, the four- 
velocity is no longer well-defined by ( 12.91 1, and the energy density w given by 
(12.12b becomes zero or negative. In order to interpret the solution, one must return 
to the Euler equations (12.6b and realize that, when w < 0, it is not possible to nor- 
malize the velocity vector, but we can still express the energy-momentum tensor 
( 12.5b directly in terms of the potential \ff, that is, 

T a p = 1VaVp -(vVyka/3- (2-18) 

Importantly, this expression is well-defined and regular for all values of y/ a \ff a ■ 
Moreover, if \ff v = \ff_ v vanishes so that w vanishes, then the only non-vanishing 
component T a p is 

T uu = 2(xi/ U ) 2 , 

which coincides with the energy-momentum tensor of so-called null dust matter. 

The regime w < is most easily understood in terms of comoving coordinates. 
Suppose first that w = y a \ff a >0so that \ff ,a is a timelike vector (due to our choice 
(+, — ,— , — ) for the signature). The comoving coordinates (T,X,Y,Z) are defined 
as follows. Set T = \ff(u, v),Y = y, Z = z, and define the function X =X(u, v) via 

dX =e b{ "- v) {xi/ v dv~xi/ u du). 

The integrability condition for such a solution to exist is precisely ( 12.141) . Noting 
that dT = \ff u du + \jf v dv, we obtain the spacetime metric in comoving coordinates 

g= (dT 2 -e' 2b dX 2 ) ~e 2b (e 2c dY 2 + e- 2c dZ 2 ) . (2.19) 

When \jf a Xjfa > 0, the variables T,X are timelike and spacelike coordinates, re- 
spectively. 

Within the above setting, we can now consider the regime \ff a Xffa < 0, for 
which T is now a spacelike coordinate and X a timelike one. Let e T ai e Y a and e z a 
be the corresponding unit vector fields so that 

XX T T Y Y Z Z 

8 = e a e p - e a - e a - e a e p . 
Then, the energy-momentum tensor ( 12.181) takes the form 

which corresponds to a matter with a positive energy density w := —iff? iffy > 
and an anisotropic stress tensor with eigenvalues w, —w, and —w. This tensor does 
satisfy the weak energy condition and, therefore, should be regarded as "physical". 

The above property can also be established by observing that, in the regime 
under consideration above, the fluid equations reduce to the one of a scalar field. 
We refer the reader to Christodoulou Hill for a detailled discussion of equations of 
state for fluids. From now on, we regard \ff as the main fluid unknown, and (12.181) 
as the main expression of the energy-momentum tensor. 
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3 The characteristic problem for metrics with weak regularity 

3.1 Normalization and choice of coordinates 

We can take advantage of the coordinate freedom to simplify radically the set of 
nonlinear equations ( 12.14b and ( 12.15b derived in the previous section. Namely, as 
we show now, the metric coefficient c and the velocity potential iff can be regarded 
as the essential variables and are governed by singular wave equations. 

Solving the third equation in ( 12.151) , that is, b m , + 2b u b v = 0, is straightforward, 
since it is equivalent to the wave equation (e 2b ) uv = 0. Hence, there must exist 
functions /, g such that 

e 2b =f(u)+g(v)>0. 
Since the transformations mi— >■ U(u) and vh>V(v) do not change the form of the 
metric ( 12.1b . we may choose the coordinates u, v to coincide with the functions / 
and g, respectively. For definiteness, we consider the case that both / and g are 
decreasing, so that b decreases toward the future and a singularity is expected in 
finite time in the future direction. It is convenient then to adopt the normalization 

/(«) = -£ ii, g{v) = -\v, (3.1) 

in order to easily recover certain particular solutions available in the literature. 
Then, the metric coefficient b is simply 

b=\\og(\\u + v\), (3.2) 

and the set \u + v < 0} is the region of physical interest, while the hypersurface 
u + v = corresponds to a (physical or coordinate) singularity. Hence, we are 
treating here the situation that b is decreasing toward the future. 

With this choice of coordinates, the spacetime (^#, g) is described by the fol- 
lowing remaining equations: 

2(m + v) yv + V« + V^ = 0, 
2(u + v)c m , + c u + c v = 0, 

au = {c 2 u + Kyl) (m + v)- ±(w + v)-\ (3.3) 



a v = (c4 + K\^){u + v)-l(u 



Wuv + ^ + frr + ,.. , 2 = 0, (3-4) 



a«v = -c u c v - Jcy«yv + ?(" + v) 2 - 

Thus, each of the functions \ff and c satisfies the same singular wave equation, 
which is a special case of the Euler-Poisson-Darboux (EPD) equations 

{u + vf 

with given constants a,j3,y. Such equations were first discussed systematically 
by Darboux [14] (Chap. Ill and IV); for more recent material, see l4l l35ll . 

Once c and iff are determined by the first two equations in ( 13.31) . one determines 
the coefficient a from the third and fourth equations. Note that the compatibility 
condition {a u ) v = {a v ) u is then automatically satisfied by virtue of the first two 
equations in ( 13.3b . The last equation for a uv in (13.3b is redundant. 

Once a,c, iff are determined, the mass-energy density w is recovered by Ber- 
noulli's law ( 12.13b . and when w > the fluid velocity is given by ( 12.9b . 



10 



Philippe G. LeFloch, John M. Stewart 



3.2 The essential field equation in a characteristic rectangle 

We now begin our analysis of the first equation (for instance) in ( 13.3b . i.e. 



which (with some abuse of notation) we refer to as the essential field equation. It 
is appropriate to pose the characteristic initial value problem. Given uq < u and 
vn < v, let P be the two-plane (u,v) and S be the two-plane (wq,vo). The plane P 
is assumed to lie in the chronological future of S, that is, u > uq and v > vq. 

Based on the past of P and the future of S, it is natural to introduce the two- 
planes R = (uo,v) and Q = (»,vq) and the associated region 



defined as the diamond-shaped region with boundary PQSRP. The value of I// is 
then specified on each of the initial hypersurfaces 



and we aim at deriving an explicit representation of y/(P) = 1/a(w,v) in terms of 
these data. Since the coefficients of ( 13.5b are singular on the hypersurface u + v = 
0, we assume that uq + vn < u + v < 0, that is £> lies to the past of this hypersurface. 
Later, we will examine the behavior of the solutions as u + v — )■ 0— . 

Remark 1 1 . Throughout this paper, one could assume < an + vo < u + v with 
virtually no change in the forthcoming analysis, and this would indeed be relevant 
for an analysis in the past direction. 

2. It is worth keeping in mind that not all solutions of (13.5b become singular on 
the line u + v = 0, and an obvious counter-example is v) = 3u 2 — 2uv + 3v 2 . 
However, iff does become singular for generic initial data, as we will show later 
(Sectionl5^i. 



3.3 Representation formula 

Our analysis of the equation ( 13.5b is based on the Riemann function approach. 
From now on, we regard the coordinates (w,v) of P as fixed and we introduce a 
new independent variable (u',v') G @. The operator adjoint to L (defined in ( 13.51 )) 
applies to functions (p = (p(u',v') and reads 



L[y]:=Yw + U u + v ) (Vh + Vv) = ° 



(3.5) 



^ = ^(mo, v o;" i v) :=SR, 



jY — c/K(an,vn;M,v) := SQ, 



L*[(p] := <p„ v - ±(u' + v')- l (<p uf + (ps) + {u' + v')- 2 f, 



(3.6) 



which is also of Euler-Poisson-Darboux type. By setting 




(3.7) 



the adjoint equation L* [(p] = is found to be equivalent to 

mV + !(«'+ V')- 1 (0„' + V ')=O: 



(3.8) 
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which coincides with the original operator ( 13.51) . 

We are going to construct a special solution of L* [9] = which satisfies a 
backward characteristic initial value problem, now with data posed on PQ and 
PR. Specifically, we choose these data to be 

/ , / x 1/2 

<p(u,v')=l — — on PQ, v < v' < v, (3.9) 



and 



/ 1 1 \ V 2 

<p(u',v) = f : - J on PR, m (3.10) 

The reason for this choice will become clear shortly. This solution depends on 
the independent variable (u',v'), as well as on the fixed parameter (w,v), and it 
is convenient to write (p = <p{u' ,v';u,v) to indicate this dependence. The solution 
(p of d3-7l) - d3T8b satisfying d 3 . 9 b - < f3TT0b is commonly called the Riemann function. 
Such a solution (p does exist since (13.8b is a linear partial differential equation 
with regular coefficients since the diamond 3>, by assumption, does not intersect 
the singularity. 

In view of the definitions of L and L* , we have 

= 2(( P L[ ¥ }- V L*[( P }) 

= ((fV)u' - 2 V<P«' + ("' + v')~ 1 ¥<P) , + ((fV)v' - 2^<Pv' + ("' + v')' 1 Wf) , 

which we now integrate over @ by using Stokes' theorem to convert the right-hand 
side to line integrals. After a straightforward integration by parts and dividing by 
a factor 2, we obtain 

(y<p)(P) - (V<P)(G) - (y<p)(r) + 

' f Q + n w{w+ v 'y\-^)dv' 

J r + J q JW (5 ("' + v'Y 1 9 - %') du' = 0. 

In view of the initial data d3.9l l assumed by (p, we see that the contribution from 
the segment QP vanishes, while similarly ( 13.10b implies that the integral over RP 
vanishes. We also use that these two pieces of initial data for (p imply (p{P) — 1. 
Thus, we find the following formula for the general solution of (13.5b 

Y(u,v) —(p(u,vo;u,v)Y(u,vo) + <p(uo,v;u,v)\if(uQ,v) 

f" 

- (p(uo,v ;u,v)y(uo,vo) - J y/(u ,v )A[(p](u ,v ;u,v)du ^ ^ 

- f \p-(u ,v')A[(p}(u ,v';u,v}dv', 

JVQ 

which is the promised representation in terms of data on J/_ and jV , with 

A[(p](u' ,v Q ;u,v) :=(p u ,(u' ,v ;u,v) - \{u' +vo)~ l (p{u' ,v ;u,v), 
A[(p}(uo,v';u,v] :=(p v i(uo,v';u,v) - ^(uo + v')~ l (p(uo,v';u,v). 
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Sufficient regularity is assumed for the time being. Later on, from the regu- 
larity of the Riemann function it will follow that the above formula makes sense 
as long as the data y/(t<o, ) an d Vt"; v o) are locally integrable. Furthermore, pro- 
vided the initial data are more regular and admit locally integrable, first-order 
derivatives, an integration by parts gives an alternative (and somewhat simpler) 
representation: 



f" 

i/a(w,v) = (p(m ,vo;m,v)v/-(mo,vo)+ / q>(u ,vo;u,v)B[\j/](u ,vo)du 

(p(u ,v ;u,v) B[\ir}(u ,v'}dv' , 



I'O 



where ^("o, vo) and 

B[v](«',v ) :=%'(«', v ) + \{u +v )~ V(«'.v ), 
B[y\ ("0, A := W (mo, v') + \ (uq + v')~ 1 \j/(u , V) , 

are given by the prescribed characteristic data. 



(3.13) 



(3.14) 



3.4 Riemann function for the essential field equation 

The formula ( 13.13b would have no utility, unless we can construct the Riemann 
function explicitly, which we now do. The results in this subsection are given in 
sufficient detail that they can be verified by (tedious) calculations. The reasons 
why the approach below works require further knowledge of the theory of the 
Euler-Poisson-Darboux equation, for which we refer to 1141 : see also [4 35] and 
the references therein. 

We start from the equation (13.8b for the "dual" Riemann function 9 and seek 
a homogeneous solution depending essentially upon the ratio u! jv'\ 

e{u',v') = (y')- l,2 y(-u'/v'). (3.15) 

Some elementary manipulation shows that 9 satisfies ( 13.8b if and only if y = y(z) 
satisfies 

z(l-z)y"(z) + (1 - 2z)y'(z) - \y{z) = 0. (3.16) 

Interestingly, this is a particular case of the hype rgeome trie equation, studied for 
instance in [1] (Section 15.5) and 1291 . 

In order to impose the boundary conditions needed to identify the function y 
we need the following result (verifiable by brute force0): if 6(u' ',v') satisfies ( 13.8b . 
then so does 

~, , ,. 1 / u' — u v' + h\ 

9(M ' V) = { ui + V )^ { V-u)^ e {-—v— V )' (3 - 17) 

2 Alternatively, one can o bserv e that {33} is the cylindrically symmetric wave equation in 
(3+1) dimensions and that d3.17t is a conformal transformation and, in fact, an inversion. 
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where, as before, u and v are regarded as parameters. Thus, taking into account 
the transformations ( 13.7b . ( 13.151) . and ( 13.171) and recalling the boundary conditions 
CLE) and flTTOl we find 

<p(u',v';u,v)= !fm ( V ,\ v, 2 y®> < 3 " 18 ) 

where 

z = zW;u, V ):= ^l^frA (3-19) 
{V + u){u + v) 

and y is a solution of d3.16b satisfying 

y(0) = l. 

Now the hypergeometric equation J3.16I) has regular singular points at z = 0, 
Z = 1, and z = °°. Near z = 0, there exist two independent solutions with asymp- 
totic forms yi(z) ~ z° and ^2(2) ~ z°logz, respectively. It is clear that the initial 
condition y(0) = 1 picks out unambiguously y\ (z), i.e., 

y{z)=F{\,\;\-z)i (3.20) 

where, by definition, F(a,b;c;z) is the standard hypergeometric function and is 
regular for |z| < 1. 

Thus, we deduce from ( 13.181) the following expression of the Riemann func- 
tion: 

<K"V;«,v) = (-^f) V2 (±±$) Uz), (3.21) 

\u + v J \u + v J L L 

Now on PQ and PR we clearly have z = and, in particular, z(Q) = z(/?) = 0. 
Along QS, the variable z increases monotonically to its maximum value z(S). Note 
that 

(u' + v')(u + v) rn 
1-Z= / , , w , , — r£ 0,1), (3.22) 
(«' + v) (v' + It) 

so that, as long as u + v remains bounded away from zero, i.e., as long as P does 
not approach the singular line u + v = 0. The condition z is obvious from 
(13.19b . while z < 1 is obvious from ( 13.22b : then z 6 [0, 1) remains bounded away 
from 1, and the Riemann function is regular. This is true, in particular, on the lines 
QS and RS. 

Inserting ( 13.21b into ( 13.13b . we arrive at the following main conclusion. 

Proposition 1 The solution of the characteristic initial value problem associated 
with the essential field equation ( 13.51) is given by the general formula d3.13b - d3.14l) 
in terms of boundary data \jf(-,Vo) and \ff(uQ,-), in which the Riemann function 
reads 

, 1 1 , /(/+v'\ 1/2 /«'+v'V /2 (\ 1 (v'-v)(m'-m)\ 

(p(u',v';u,v)=(- FlLl',1',)- J ±- '-) (3.23) 

rv ; \u' + vj \u + v' J V (v' + u)(u'+v)J 

and the (hypergeometric) function F is singular when its last argument tends to 1. 
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3.5 Definition and existence of W ' regular spacetimes 

We now return to the full set of Einstein equations (13.3 b or, more precisely, 

2( M + v)i// MV + y„ + i// v = 0, 

2(u + v)c uv + c u + c v = 0, 

a u = (c 2 u + k yl) (u + v) - i (u + v) " 1 , ' " ' 24 ' 

a v = (cv + k {u + v) - \ (u + v) " 1 , 

which we refer to as the reduced Einstein equations for plane symmetric space- 
times (y$,g) (in the chosen gauge). Now, we are interested in encompassing solu- 
tions lff,c,a with weak regularity. Specifically, we propose to search for solutions 
\j/,c in the Sobolev space W 1 ' 2 of functions which are square-integrable together 
with their first-order derivatives. This class is natural since the curvature is then 
well-defined in the distributional sense, and the Einstein equations ( 12.2b hold as 
equalities between distributions in the dual space W~ 1 ' 2 , as established in H231 . 

Given any uo < u and vo < v with u + v < 0, we focus on the characteristic 
initial value problem within the (regular) region 3 = S>(uo, vo; u, v) C ./# defined 
in Section 13.21 The boundaries jV_ and JV are null hypersurfaces on which we 
prescribe data and which intersect on the two-plane 

,9> = ^*(m ,vo) :={u = u ,v = vo}. 

Without loss of generality, we can normalize the metric coefficient a to satisfy 

a(wo,vo) =0 on the two-plane 3 s . (3.25) 

At this stage, it is most convenient to provide a first description of our results in 
terms of the (coordinate dependent) coefficients a,c and potential \ff. But, later in 
Section |5](cf. Theorem[2]l after some further analysis, we will restate these results 
in a fully geometric (i.e., chart invariant) form. 

It is convenient to introduce the following notation for any function / = f(u, v) 



MZ, V0 \f}-= ^p ([ U \f(;v')\ 2 du'), 
f]:= sup ( r\f(u',-)\ 2 dv') l/2 



T7 U - V 



and, by extension, 

M£XLfl == (Jjf(u',v )\ 2 du') 1/2 , MT Q ; Q [f] := (£ \f(u , v')\ 2 dv') "\ 
Other straightforward extensions of this notation will be used, for instance 

MZ V0 [fA=m, V0 [f}+M%:, fs [h] 



if two functions /, h are given. 
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Definition 1 (Notion of W 1 2 regular spacetime) Given uq , vn and u , v with u < uq , 

v < vo and uq + vn < 0, a W 1 ' 2 regular spacetime satisfying the Einstein equations 
d3 .24b in the characteristic rectangle 3 = &(uq , vo; u, v) and describing self-gravi- 
tating irrotational fluids, is determined in characteristic coordinates (m,v) by three 
continuous metric coefficients a,b,c (cf. J2.ll) ) and a continuous fluid potential iff 
(cf. (E3 and d2TT2l) ) such that: 

1. The coefficient is given by the explicit formula ( 13.2b . 

2. The derivatives of a,c, l// are differentiable in a weak sense and the (semi-) 
norms 

M%[ V0 [\a u \ l/ \c u , ¥u ], NC;°[\a v \ x l\c v , Wv ] 

are finite. 

3. The reduced field equations ( 13.24b hold in the sense of distributions. 

4. The function a satisfies the normalization ( 13.25b . 

Our terminology "W 1 ' regular spacetime" is motivated by the fact that c and 
\ff are the essential unknowns of the problem and, by our definition, have square - 
integrable first-order derivatives. As we observed earlier, initial data are required 
for the functions c and iff, only, while the function a can be recovered afterwards. 

Recall that the potential iff determines, both, the velocity field u and the mass 
density w. In view of ( 12.131) , the regularity assumed in the above definition implies 
that the spacetime integral of the mass density 

// \w(u',v')\du'dv' 

J 'J V(ho,v ;h,v) (3.26) 
< Ae~ 2 ^ a (« - a) 1 / 2 (v - v) 1 ' 2 m^WX ,«WA 

is also finite. On the other hand, no general estimate is available for the velocity 
field, since the norm of the gradient of \ff, that is, w, may well vanish, at which 
point ( 12.91 ) is ill-defined. Recall that iff is the main fluid variable, and d2.18b is the 
primary expression of the energy-momentum tensor. 

Theorem 1 (Well-posedness theory) Let 2> = S)(uo,vo;u,v) be a characteristic 
rectangle that does not intersect the singularity hypersurface. Let iff, c and 1/7, c be 
continuous functions (of a single variable) defined on the null hypersurfaces 

jV_ = {«o < u' < u; v' = vo}, ^V={u'=uq; Vo <«'<«}, 

respectively, and satisfying the continuity conditions ^("o) = V^(vo) and c(uq) — 
c(vo), such that the semi-norms 

MSo^vg, <;; [cv,¥v] 

are finite. Then, there exists a unique W ' regular spacetime (in the sense of Def- 
inition^ determined by functions a,b,c,lff : S> — > K that satisfy the reduced Ein- 
stein equations ( 13.24b of self-gravitating irrotational fluids and assume the initial 
data 

c ('jVo) = c, y/(-,vo) = \ff on the hypersurface JV \ 
c(uo,-) = c, Xff(uQ,-) = \j/ on the hypersurface jY '. 
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Some additional remarks about the regularity of our solutions are in order. 
Recall the embedding W 1 - 2 C C 1 / 2 valid on one dimension, so that any solution is 
Holder continuous in each characteristic variable. The second-order derivatives of 
the solutions iff, c (which can always be defined in the sense of distributions) need 
not be functions in a classical sense. However, the mixed derivatives Xff uv and c HV 
do have some regularity and, specifically, the spacetime integral 



9>(u Q ,v Q \u,v) 



is also finite. 



|V«v| 2 ) du'dv' 



(3.27) 



Proof. Given the reduction analysis and the observations already made in the 
present section, the proof is now direct. Indeed, we have derived the integral for- 
mula (I3.13b - d3.14b . which provides an explicit expression for the solutions c to 
the first two equations in ( 13.24b in term of their characteristic data \ff,c and 1/7, c 

prescribed on the null hypersurfaces J/_ and jV , respectively. Since the Riemann 
function is bounded (at least) in the domain under consideration, as explained 
before Proposition Q] all the integrals in ( 13.131 ) make sense since B[\ff] and B[\jf\ 
are (in L 2 and thus) integrable. The existence of functions \ff, c : 21 — > K satisfy- 
ing the reduced Einstein equations in the distribution sense and the desired initial 
conditions is now clear. 

Furthermore, estimates on the norms of these solutions, as required in Defini- 
tion [TJ can be established from the integral formulation, as follows. Considering 
for instance the function for all («',v') G 2(uq,vq;u,v) we have 

\lf(u',v') =(p(u',vo;u',v') \jf(u') + (p(uo,v;u' ,v')\j/(y') 

- <p(mo,vo;m',v')3 (V("o) +F( v o)) 

- f" w{u")A[(p]{u",v Q ;u'y)du" - f xj/(v")A[(p](uo,v";u' ,v')dv" . 

Using that the Riemann function is regular away from the singular hypersurface, 
we compute the derivative \\f u i and obtain 

I v«'("'/)l < I v^OOl + + Iv 7 ^')! + \v( u o)\ + l¥(vo)| 

+ ( \Y{u")\du"+ f \xj/(v")\dv", 

in which the notation / < g means that there exists a constant C > such that 
/ ^ Cg. Since a direct calculation allows us to control the sup norm of \ff within 
the spacetime region, we obtain a pointwise control of the u' -derivative of iff in 
terms of the same derivative of the initial data: 



\Yu'(u',v')\< |V^,0')I+ sup + 

3/{itQ,v a ;uy) 
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By integrating this inequality along an arbitrary characteristic line allows us to 
control the desired W 1,2 norm of the solution with the same norm of the initial 
data: 

mz vo ivu}= sup ([ u \ ¥u (.,v')\ 2 du') 1/2 

<M^°v fc]+ (\y\ + \W\)- 

3>(u ,v ;u,v) 

The same arguments apply to \ff v , as well as to the coefficient c. 

Next, we return to the key set of equations ( 13.241) and determine the metric 
function a as follows. By integrating the relevant equations in ( 13.241) along the 
initial hypersurfaces Jf_ and jV , respectively, and using the normalization ( 13.251) 
we obtain the initial values for the function a: 

a(u') := a(u', v ) = f ((c? )M + k£)(u") («" + v„) - ^(«" + v )- 1 ) du" , 

a(V) :=a(uo,v')= [ {{c\ v + ^)(v") («" + v") - \{u" + vV) dv" . 

Both expressions provide bounded functions since the prescribed characteristic 
data belong to W 12 . Having determined a on the two null hypersurfaces SR and 
SQ, we can use once more the last two equations in ( 13.31) and, by integrating along 
characteristics, determine the function a within the whole spacetime region 2$. The 
equations ( 13.24b show that the functions a u and a v have the same regularity as the 
products c u \if u and c v lff v , respectively, that is, are integrable along characteristic 
lines. Hence, M^v [Kl 1/2 ] and M^°[\a v \ l l 2 \ are finite. □ 



4 The characteristic problem for NP scalars with weak regularity 

4.1 Choice of the tetrad 

Our objective here is to reformulate more geometrically the characteristic initial 
value problem when data are given on two intersecting null hypersurfaces and, 
specifically, to identify which physical data should be imposed on these hypersur- 
faces. To emphasize that the null coordinates under consideration can differ from 
the ones constructed in the previous section, we denote them by (U, V). We are 
going to re-derive the expressions of the field equations (obtained in Section[2]) via 
the formalism introduced by Newman and Penrose 0281 . following here the nota- 
tion in 1 34 1 . Our main conclusion in the present section, as stated in Proposition [2] 
below, is a restatement — in the NP notation — of the weak regularity conditions 
introduced in Definition [T] 

The tangent space can be described by a basis n a , l a , m a , m a of complex- 
valued null vectors satisfying the normalization: 

1 = l a n a = n a l a = -m a m a = -m a m a , 
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while all other contractions vanish. In view of the expression (ED of the metric 
for plane symmetric spacetimes, we choose the tetrad 




in which the treatment in the U - and V-directions is not symmetric as this has 
some advantages in the applications (to conveniently prescribe the incoming radi- 
ation). As is standard, we use the notation D, A , 5, 5 for the directional derivatives 
associated with the null tetrad. 

We begin by describing general formulas satisfied by Newman-Penrose scalars 
in such a null frame and, next, will specify our particular choice of tetrad. 



4.2 Expressions for the Newman-Penrose scalars 

In the NP notation, the Christoffel symbols associated with the metric (12.1b are 
represented by the following twelve connection NP scalars: 

a = P = y = 0, £ = au, 
K=0, X=e 2a c v , 

; (4.i) 

H=e by, v = n = 0, 

p = -bu, a = —cu, T = 0. 

while the Ricci curvature is represented by the following seven Ricci curvature 
NP scalars: 

#oo = e~ Aa (buu ~ 2aubu + b% + c%) , 
*bi = #12 = 0, 

#02 = e~ 2a (c uv + bucv + b v cu) , 

(4.2) 

i> n = \e~ 2a {-a uv +b v b v -c u c v ) 7 

#22 = e Aa (-b vv +2a v b v -b\-c\), 

A = \e~ 2a (auv + 2buv + 3bub v + cuc v ) , 

and the Weyl curvature (or free radiation) is represented by the following five Weyl 
curvature NP scalars: 

% = -cuv + 2(au-bu)cu, 
W l= W 3 =0, 

^2 = 3 e~ 2 " {-aw + buv +2 Cu c v ) , ^ ' 

% = e Aa (-cvv + 2{a v - b v )c v ) . 

Of course, the Ricci curvature is related to the matter tensor via the field equa- 
tions 

#aj3 = -f (jap ~\T gafi), (4 4) 

A _ K rp rp rji(X 

yi — 24 ' * — « * 
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To make this more explicit, it is convenient to decompose the fluid velocity vector 
u a in the form 

u a =Zl a + Wn a , Z = u a n a 7 W = u a l a , (4.5) 

so that the condition u a u a = 1 becomes WZ = j. From the expression of the 
energy-momentum tensor (for general pressure p) 

T a P = + p)u a uP -pg aP , T=w-3p, 
and in view of ( 14.4b we can compute 

*a/3 = -j{w + p)(u a up - \g a ft), 
A = -±(w-3p). 

We thus find 

<£o2 = 0, (4.6) 

as well as the relations 

*bo = -f (w + p)W 2 , <t> n = -|(w + p), *22 = -f (w + p)Z 2 , 
A = -±(w-3p), 

which, after imposing p = w for a null fluid, become 

<%, = -kwW 2 , <£n = -fw, <P 22 = -kwZ 2 , 
A = ±w. 

In view of the third relation in (14.21) , the equation ( 14.6b is equivalent to a wave 
equation for the function c 

cuv +bucv +b v cu =0, 
which allows us to recover the key equation ( 12.16b derived earlier. 



4.3 Evolution equations for the NP scalars 

Before we can proceed further, we need to identify which variables are the essen- 
tial dependent variables that require initial data. First of all, it follows from (14.7b 
that <i>oo , <i>22 , 3*1 1 are not independent but satisfy the constraint 

Ci := (*oo 3* 22 ) 1/2 +2 4*! i =0. (4.8) 

(Recall that 4>oo and 4>22 are both real and have the same sign.) We thus impose the 
constraint C\ = everywhere and regard <E>n as a redundant dependent variable. 
For if we know <£>oo and 4>22 at a point where w ^ 0, then (4.8) determines 4>i 1 up 
to a sign, and (4.7) requires <Pu ^ 0. 

Similarly, <£>ii and A are not independent but are related algebraically: 



C 2 :=3*n+3A =0, 



(4.9) 
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and thus A can be viewed as a redundant dependent variable. Interestingly, the 
relation ( 14.9b is equivalent to an evolution equation for b, i.e. 

bvv+2bub v = 0, 

which allows us to recover the equation (12.15b . 

Third, the NP scalars satisfy one further algebraic constraint because dp = 
8g = 0. By taking into account the property C\ = C2 = this constraint can be 
written in the form 

C 3 :=V 2 + 2A-(pn-o?i) =0, (4.10) 

which means that W 2 can also be regarded as a redundant dependent variable. 

We are now in a position to list the complete set of evolution equations for the 
following nine non-redundant NP scalars 

<£()0, *22, %, %■ 

First of all, the field equations consist of four evolution equations in the D- 
direction: 



D 3 
D 4 



= Dp-(p(p+2£) + CT 2 + <£ o) =0, 

= £>CT-(2(p + £)CT + fo)=0, 

= DX - ((p - 2e)A + ap) = 0, 
= Dii-2{p-e)!X = 0, 



(4.11) 



and five evolution equations in the A-direction: 



A 


= Ae- 


- (-p/i + oX+6A) = 0, 


Ax 


= Apn 


-2pp =0, 


A 2 


= Ao- 


\-jio + ?ip =0, 


A3 


= Ak- 


F2juA + % = 0, 


A 4 


= Ail - 


hjU 2 + A 2 + <J> 22 = 0. 



(4.12) 



Note that the constraints C\ = C 2 = C3 = have been imposed in the equations 
above, and will be imposed as well in the Bianchi identities below. The relevant 
two contracted Bianchi identities 



D 5 :=D<£>22+4 J u'£ii +2(-p + 2£)<£ 2 2 =0, 
A 5 :=A<£oo + 2jU<&oo-4p<i>ii =0, 

while the two independent Bianchi identities reduce to 

D 6 := DW 4 - (-3Af 2 + (p - 4e)f 4 - 2A4>n + o<P 22 ) = 0, 
A 6 := AW - {-LiVo + 3oV 2 - A4> 00 + 2ff*n) = 0. 



(4.13) 



(4.14) 



It is important to notice that the quantities D n and A„ (n = 1 , . . . , 4) are not 
independent but satisfy four differential identities. To obtain these, observe that 
for any scalar / 

D(Af)-A(Df) = -2eAf, 
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showing that D and A do not commute in general. We can for instance express this 
general identity for the scalar p, and then rewrite Dp in terms of D\ and Zip in 
terms of A\ . Then we replace all D-derivatives of NP quantities by the appropriate 
and all A -derivatives by the corresponding A^. Finally, making use of the three 
constraint terms Q we arrive at the desired identity and, applying the same process 
to a, X and /i in turn, we obtain 

DAi-ADi =2pA + 2fiDi + 2pA\ +2oA 2 + 2pD 4 + A 5 +4pC 2 , 
DA 2 - AD 2 = 2a A + XD\ + 2oA x + \iD 2 + 2pA 2 + pD 3 + oD 4 

+ A 6 + 2aC 2 + 3aC 3 , 
DA 3 -AD 3 = -2?iAQ+?iAi++iiA 2 + 2ixD ? , + (p-4e)A 3 +2?iD4 

+ (744 + £>6-2AC2-3AC 3 , 
DA 4 - AD 4 = -2ju z\ + 2\iA x + 2XD 3 + 2\iD 4 + 2pA 4 +D 5 - 4/xC 2 . 

Finally, once the NP scalars have been determined, the metric coefficients in 
(12.1b are simply recovered by integrating out the following equations for a,b,c: 

au = e, 

bu = ~p, b v =e- 2a pL, (4.16) 
cu = — a, c\ = e~ la X. 

Observe that there is no equation for ay. (Recall that the coordinates u and v do 
not play symmetric role, due to our choice of tetrad.) 



4.4 Revisiting the characteristic initial value problem with regular data 

We are ready to formulate the characteristic initial value problem for plane sym- 
metric spacetimes. We are given a point (Uq, Vo) representing the two-plane 

9 := {U = U , V = V }, 

and consider the two hypersurfaces J/_ = {V = Vq} and JV = {U = Uq} with 
intersection . 

The null coordinates and the tetrad are geometrically defined as follows. We 
first choose the vectors / and « at ^ to be future-directed null vectors normalized 
by the condition g([,n) = 1 and tangent to jV_ and jY , respectively. The vectors 
are determined up to the transformation (/,«) i-> (p[,n/ p) for p > 0. The quotient 
manifold Q of the four-dimensional spacetime by the symmetry group (i.e. the 
two-dimensional Euclidian group) is represented by each of the timelike surfaces 
which are orthogonal to the group orbits. The vectors /,« can then be extended 
as geodesic fields in the future direction and, in the quotient picture, their integral 
curves are nothing but ^K, jV, respectively. We then define the coordinate U on 
jY to be the affine parameter of / normalized so that U = Uq at the coordinates 
V is defined similarly from n. Next, the functions U,V can be uniquely extended 
to Q by the requirement that their level sets are outgoing and incoming null curves, 
respectively. 
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From the gradient of these functions, we can then determine l a = g a ^djjV and 
g al3 dpU, which are null geodesic fields satisfying the orthogonality con- 
ditions l a d a V = = n a d a U . Hence, I is tangential to jV_ and n is tangential to 
JV . Consequently, along J/_ we can write l = Cl for some constant C > 0, while 
along JV we have n = Cn for some constant C > 0. Then, we write g(l,n) = 
g a Pd a VdpU = g uv = e~ 2a , and we define V = e 2a l and ri = e 2a n. It follows that 
l' a d a U = g(l',n) = e 2a g(l,n) = 1 and, similarly, n' a d a V = g(nf,l) = e 2a g(n,l) = 
1. Consequently, we actually have V = dy and «' = cfy, and it follows also that 
l a d a U = e~ 2a and n a d a V = e~ 2a . 

Now, along the hypersurface jV_ and in view of l a d a U = e~ 2a and l a d a ll = 1, 
we conclude that C = e 2a along Jif . An analogous argument tells us also that 
C = e 2a rJV and, by continuity of the function a at it follows that for some 
k > we have C = C = k and therefore e~ 2a =k along j¥\} jV . Recalling that 
/ = C7, n — Cn, and C = C = k, we get Z = (1 /k)l and « = (l/k)n, so that 

1 =g{Ln)=k- 2 g{l,n)=k~ 2 e~ 2a = l/k. 

We conclude that /c = 1, which justifies the normalization that a = on JT}J>jV . 
Our formulation is based on prescribing the following NP scalars 

p,a,A,jU on the planed, (4.17) 

together with data 

a , ff) , <J>oo on the hypersurface ^/K , (4.18) 

and 

a, % , <P22 on the hypersurface jV . (4. 19) 

These data can be given freely, except that we require <J>oo ^ on jV_ and <J > 22 ^ 
on JV , in agreement with J4.7I) . In addition, since regularity is required in the 
present section, we impose that 

a(Uo)=a{V ). 

Furthermore, without loss of generality we impose the normalization 

a = b = c = on the plane 9 s . (4.20) 

Our aim is to determine the solution components 

e,p,a,X,ii, *oo,*22,f / o,f4 

in the future @ of JV_ U ~ZV . 

Consider first the solution on the two-plane . In view of (14.17b and ( 14.181) 
we know £, p, a, X and ji. Also since we know <Pqq and <£>22 at ( 14.8b fixes 
^>ii ^ 0. The scalar A follows from ( 14.9b , and W2 from J4.10I) . This determines the 
complete set of dependent variables on 9. 

Next consider the solution on the hypersurface jV . Since ci(-,Vq) — a is pre- 
scribed, we know £ from ( 14.16b . Next, ( 14.1 lb gives a coupled system of nonlinear 
ordinary differential equations of Riccati type for (p , a), with known source terms 
and initial data known on the plane & > . Hence, using the prescribed data at (Uq , Vq) 
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we see that this system has a unique solution (p,o) defined on a maximal inter- 
val, denoted by [Uo,Uq), Returning to ( 14.1 lb . we may then solve a linear system 
for the unknowns (A,ju) on l A_ with prescribed initial data on 2? . Next, using the 
condition C\ = to express <P\\ ^ in terms of the known <Pqq and the unknown 
4>22 we can regard D5 = as a (sub-linear) ordinary differential equation in U 
for <E>22; hence, solving this equation delivers <E>22 and hence <Pn and A. Then, 
the condition C3 = produces W2. Finally, after substituting for W% and <Pn, the 
condition = is a linear equation for W4 with known data, and so we have 
determined a complete set of dependent variables on jV_. 

When we consider the solution on the hypersurface jV a similar argument 
applies. We see that A\ = A2 = At, = A4 = is a coupled set of Riccati ordinary 
differential equations for X , /I with known sources and initial data on 3 s . Using the 
data at (Uo, Vb), we see that this differential system has a unique solution defined 
on a maximal interval [Vb, Vo)- We then solve a linear system in p, a. Using C\ = 
to express 4>\\ ^ in terms of the known 4>22 an d the unknown c&oo we see that 
A$ — is a (sub-linear) equation for <Pqo with initial data on 8?, This delivers <i>oo 
and hence 4>n and A on jV for a suitable finite V-interval. Then C3 = produces 
W2. Next Af, = is a linear equation for Wq with known source terms and given 
data at 3, which delivers Wq. Finally Aq = is a linear equation with known initial 
data for £ and hence we have determined a complete set of dependent variables on 
Jf for a suitable finite V-interval. 

The construction so far has produced solutions for £, p, a, X, jU, <Pqq, 4>22, 
Wq and W4 satisfying D n = on .jY_ and A n = on Jf . We have available 6 in- 
equations and 7 A -equations to determine the solution for these nine dependent 
variables in the future of /U/. We note that both D\ = and zli = could 
be used to evolve p and so we choose to discard one of them. As we shall see, it 
does not matter which. We do the same for a, A and fi dropping either D n = or 
A„ = for n = 2, 3,4. We now have 9 equations for 9 dependent variables. We are 
not allowed to choose a(u,v) freely in the interior, although we know its value on 
jV . However one of the equations (14.161) implies Do = where 

D Q :=Da-e, (4.21) 

which we adjoin to our set, so that we have 10 equations for 10 independent vari- 
ables a,e,p,<7,X,H,<pQ(),<p22, i I'o,W4, and there is precisely one D-equation or a 
A -equation for each variable. These form a first order quasilinear system, pro- 
vided we use the constraint equations Q = to eliminate <P\\, A and W%. Further, 
writing the equations in the stated order it is obvious that the system is diagonal 
and each entry has either a positive /-component or a positive n-component, and 
this guarantees that the solution exists and is unique. 

The one remaining difficulty is that our system is not unique, because of the 
discarding process above — in fact there are 2 4 such systems. However, in each 
such system we set precisely four of the {D„,A n } (n = 1,2,3,4) to zero (one for 
each choice of n), so there is no a priori guarantee that the others will also be zero. 
We now examine the set ( 14.151) , setting to zero all of the Cj to zero, and those of 
the D„ and A„ which we know to be zero. We also replace the coefficients in this 
set by the solution we computed in the previous paragraph. We are left with a 
linear symmetric hyperbolic system for the remaining D„ and A n with zero source 
terms and trivial data on i A_ or jV as appropriate. Clearly the solution is the trivial 
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one. Thus the procedure described in this section generates a unique and complete 
set of dependent variables, at least for sufficiently regular solutions. 

4.5 Weak regularity of the NP scalars 

The material in the present section provides a method to solve the characteristic 
initial value problem and establish the existence of a solution within any charac- 
teristic rectangle 3(uq,vq;u,v) avoiding the singular line. Such a local existence 
result follows from standard theorems, in view of the symmetric hyperbolic form 
of the equations exhibited above. However, the above presentation does not pro- 
vide a global existence result. In that sense, the analysis in Section [331 was more 
precise and led us to an actual proof of existence, based on non-physical data 
though. Importantly, the present section has allowed us to identify the physically 
relevant data, and, in the next section, we will put together our two approaches. 

It remains to discuss the regularity of the NP scalars when the spacetime is 
solely weakly regular. By comparing the regularity in Definition [TJ with the ex- 
pressions of the Ricci and Weyl scalars derived in the present section, we arrive at 
the following regularity results. 

Proposition 2 The Ricci and Weyl NP scalars associated with the weakly regular 
plane symmetric space times defined in Theorem\l\within a characteristic rectan- 
gle 3) have the following regularity along the hypersurfaces jV_ and jY; 

<ieff u (/), w ew- h2 {JO, ^ooez^UD, (4.22) 

and 

«eiv u (/), ^er 1 ^/), <P 22 eL\^r), (4.23) 

as well as the following regularity in the spacetime: 

Moreover, these spacetimes satisfy <Pq\ = <E>02 = 3>i2 = and W\ = — 0. 

This result provides us with a geometric formulation of the regularity of W ' 
weakly regular spacetimes. 

Remark 2 In the coordinates chosen in Theorem [TJ the coefficient b is defined 
by ( 13.2b and is thus smooth, but this regularity property is tight to our choice of 
characteristic coordinates. 



5 Global causal structure 

5.1 Main result 

This section is devoted to a proof of a global and fully geometric result which 
goes well beyond the local result given earlier in Theorem [TJ We still prescribe 
data on two null hypersurfaces intersecting along a two-plane, but we are no longer 
working within a given characteristic rectangle and seek for the global structure of 
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the future development of the given initial data set. Following the discussion in the 
previous section, we consider arbitrary null coordinates, denoted below by (U,V), 
which in general differ from the coordinates (u, v) constructed in Section[3] Recall 
that, throughout, we restrict attention to plane symmetric data and spacetimes. 

Definition 2 An initial data with weak regularity consists of the following pre- 
scribed data. Let 

{JL, eHUdydz) , (JF, e"dVdydz) 

be two plane symmetric 3-manifolds (endowed with volume forms) with bound- 
aries identified along a two-plane and parametrized for some (Uq, Vo) as 

J^:={U^U }, ZV:={V^V }, &>:={U = U ,V = V }. 

1. Suppose thata,a are absolutely continuous, i.e. the integrals 

I (\a\ + \dua\)e^dU [_(\a\ + \d v a\) e"dV 

are finite, and are normalized so thata\g» =a\&> = 0. 

2. Let ^rjif^oo an d ^4^22 be (plane-symmetric) functions defined on the hy- 
persurfaces ,jY_ and Jf , respectively, with ^ <& o 6 L 1 (^£1) an & ^ ^22 6 
L l (,yV), i.e. the integrals 



I Gao&dU, [_<P 22 e a dV 

JjV J.jV 



are finite, and thatVja G W~ 1;z {jV) and W 4 e W~ x;i {^Y), i.e. = 



andW 4 = d v w[ l) with 



\^\ 2 e a -dU, 



<P { £\ 2 e a dV. 



3. Finally, one also prescribes the connection NP scalars po, Oq, Xq, on . 



Theorem 2 (Global causal structure of plane symmetric matter spacetimes) 

Consider an initial data set with weak regularity determined by (jV, e-j and 
{jV ,e a ), a plane (^,po,Ob,Ao,jUo), and prescribed Ricci and Weyl NP scalars 

(S>,*oo) ^d (V 4 ,& 22 ). 

(1) Then, there exists a unique W 1 ' 2 regular spacetime (,/df,g) determined by 
metric coefficients a,b,c and matter potential iff which is a future development 
of the initial data set satisfying the Einstein equations ( 13.3b for self-gravitating, 
irrotational fluids with the initial conditions 

(p,o,A,,n) = (p ,c7o,Ao,Mo) on &>, (5.1) 

(a, 0qo) = (^j^Oi^oo) on tne nu tt hypersurface jV , (5.2) 
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and 

(a, f^, <E>22) = (a,f4, <E>22) on the null hypersurface jV . (5.3) 

(2) The constructed development of the initial data has past boundary 

{U >U> Uq;V = V } U {U = U ;Vo > V > V } C ^U^ 7 

and, for generic initial data (in the sense defined in (3), below) the curvature 
blows up to (and makes no sense even as a distribution) as one approaches its 
future boundary 

^ Q :={F(U)+G(V)=0} 

for some functions F : G in W 1 ' 2 (i.e. having two derivatives in L l ), so that the 
spacetime is inextendible beyond S3§ within the class ofW 12 regular spacetimes . 

(3) This result holds for generic initial data, in the sense that arbitrary data 
can always be perturbed in the natural (energy-type) norm so that the perturbed 
initial data do generate a singular spacetime whose curvature blows-up on 

The coefficients e- and e a , modulo a conformal transformation, could be cho- 
sen to be identically 1 on the initial hypersurface, so that two main degrees of free- 
dom remain on each of the two initial hypersurfaces. Theorem [2] can be seen as 
a statement of Penrose strong's censorship conjecture for plane symmetric space- 
times with low regularity. The proof of this theorem requires a sufficient knowl- 
edge of the singularities of the Riemann function, which we discuss below. We 
emphasize that the theory developed in this paper applies to the matter model 
(12.18b described at the end of Section 12.11 which is equivalent to the model ( 12.5b 
as long as the energy density w = V a i//'V c( i// remains positive. 

5.2 Passage from metric data to NP scalar data 

We first discuss the "existence part" in Theorem [2] In principle, this result fol- 
lows by applying Theorem Q] and patching together local solutions constructed in 
characteristic rectangles. The main difference between the two statements lies in 
the formulation of the initial data and, thus, we need to check that the initial data 
posed on the NP scalars are sufficient to determine the initial data posed in terms 
of metric coefficients as was required earlier in TheoremQ] 

We are given initial data in terms of conformal factors and curvature NP scalars 
prescribed on the two initial hypersurfaces, as stated in J5.2b - d5.3b . To recover the 
earlier description of the characteristic data in terms of the metric coefficients, we 
first observe that the two choices of null coordinates need not coincide. So, we 
search for new characteristic variables 

u:=F{U), v:=G{V) 

for some function F, G, that remain to be identified and may depend on the pre- 
scribed data. Roughly speaking, imposing Wo and <Pqo on Jf_ is analogous to pre- 
scribing c and b, respectively and similarly, imposing W4 and <E>22 on is analo- 
gous to prescribing c and b, respectively. 

Consider first the hypersurface jV_ on which we are given 



(a,fo,<£oo) = (fl,^o,*oo). 
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Our first task is, along jV , to solve a Riccati matrix system in the variable U for the 
unknown vector (p,0_) := (p,a)|. / f. Using the quantities D\ and D2 introduced 
in ( 14.11b we find — 

Dp -p(p + 2e)-o; 2 + *oo = 0, 
Da-2(p + e)a-W J3 = 0, 
in which £ = ay = Ojj is a prescribed data along jV, thus 

<§H e+ ^i)(§) + (t) 

The solution, in general, blows-up at some finite value denoted by C/ which could 

be estimated by writing a matrix Riccati equation for P : = ^ . Importantly, 

given the regularity of the Ricci scalar <£ 00 G L l and ^ G W~ 12 we see that 
p G W 1 ' 1 G^Q and a G L 2 (^/jQ. 

Having identified the scalars p,ff on jV_ we can recover b and c by integration 
of dOT) : 

/•£/ /■{/ 

b(U) = - p(U')dU', c(U) = - a(U')dU', (5.4) 
7£/ - .7c/ 

where we have chosen the normalization b(Uo) = c(Uq) = 0. Note that b G W 2,1 [ jV ] 
and c G W ' (jjQ. We can then pursue the construction of the data on as was 
explained in Section|4]and, in particular, we recover the fluid potential as well. 

Similarly, along JV the same argument produces the value Vo and the initial 
data b,c,lff. 

The analysis made earlier to show the existence of a weak solution within a 
characteristic rectangle applies to show the existence of a solution a,b,c, \ff, where 
now the function b is not normalized a priori, but at this stage e 2h may be a general 
solution to the wave equation. 

The largest possible domain of interest is [Uo,Uo] X [Vo,Vo]- However, not 
all of it is relevant since, in general, a blow-up in the function b will take place 
before one can reach the boundary of this domain. The future boundary of the 
spacetime is determined by the function b, as follows. We determine the functions 
F, G by considering b along the initial hypersurfaces and jf . Since b is a 
geometric invariant (related to the area of the orbits of symmetry), it satisfies the 
wave equation derived earlier (in coordinates), and we can write 

e W) = _i F(c/) _i G(y) 

for some functions F, G. These functions, when increasing, are used to define a 
change of (null) coordinates, defined by 

u:=F(U), v:=G(V). 

In terms of the initial functions b and b (already computed in ( 15.41) ), we find 



F(U) = 2e 2 ^+G{V ), G(V) = 2e 2biy) +F(U ), 
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in which F(Uq),G(Vq) are arbitrarily fixed with, since b(Uo) = b(Vo) = 0, 

F(Uo) + G(y ) = -2. 

Finally, we apply Theorem[T]with suitable family of characteristic rectangles cov- 
ering the whole domain 

■= {U.Q > U > Uq;Vo >V> V q ;F{U) + G{V) < 0}. 

To complete the proof of Theorem [2] it remains to investigate the nature of the 
future boundary 

:= {F{U) + G(V) = 0} = {u + v = 0}. 



5.3 Blow-up behavior of the Riemann function 

We now return to the notation introduced in Section [3T2l within a characteristic 
rectangle, and reconsider the solutions c, iff constructed in Theorem [TJ We need to 
investigate their behavior as the "vertex" P = (u,v) of the characteristic rectangle 
approaches the singular line, i.e., u + v — > 0— . Recall that z = on PQ and PR. 
Next, consider the line QS and observe that F(j,j,1;z) = 1 at Q. In view of 
( 13.22b . when u + v — > 0— one has also z{S) — > 1— . 

The Riemann function involves the function F(j,j;l;-), which can be com- 
puted thanks to the following Euler's formula (valid for \z\ < 1) 

^hb^)=J tl/2{1 _ t y /2{1 _ zt) ^ (5-5) 

for which we refer to [29] (Chapter 5, equation (9.01)) or [ 1 1 (equation (15.3.1)). 

Clearly, as z — > 1 the denominator in ( 15.5b approaches f '/ 2 (1 — t), leading to 
an integral which diverges logarithmically at the upper end. Thus, for z close to 1, 
the dominant contribution to the integral comes from the upper end, and we can 
quantify this property as follows. By defining the function 

[t,Z) Jo (1-t)V2(i_ zt) i/2' 
we can check that, at the leading order, 

7fF(i,I;l;z)~/(l,z) asz^l-. 
On the other hand, the integral /( 1 , z) can be computed explicitly and we conclude 

F(i,i;l;z) log(l-z) asz^l-. (5.6) 

We return to the representation ( 13.131 ), and note first that in each of the inte- 
grands, the Riemann function (p has constant sign so that no cancellation can take 
place within each term. Consider, for instance, the first integral (over the segment 
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SQ) and assume that B[\ff] (• , vo) is a "generic" function so that the behavior of the 
first integral term in ( 13.13b can be determined by studying 

/( „, v) := r (•<£») " 2 "> h . ; n f'rlfr-l ) * 

y ' Ju \u' + vJ \u + v J V (v' + u)(u' + v)J 

(5.7) 

in which we have specified the relevant value of the argument z. 

In the region u + v < the function / is regular, and we are interested here in 
the limit 

u + v ■= e — > 0-, 

in which the coefficients of the Euler-Poisson-Darboux equation ( 13.51) becomes 
singular. Clearly, a singularity in f(u,v) must arise at one (or both) endpoint(s) of 
the integral. 

Near Q, u' « u, z ~ 0, and both the second and third factors in the integrand 
d5.7l) are approximately unity. Then f(u,v) picks up a term Im + vop^w + vI 1 / 2 
which is itself finite in the limit — £ = u + v — > 0. However its first derivatives are 
0{e~ x l 2 ), singular as £ — > 0. 

Near S, u' w Mo and we see from ( 13.221) that 1 — z= 0(e) and so F(|, |; l;z) ~ 
— 7T _1 log£ and f(u,v) ~ — ^Tlog£ where is a strictly positive constant. Not 
only does f(u,v) become singular when £ = w + v — )■ 0, but its first derivatives 
are 0{e~ x ). This behavior dominates the milder singularity which originates from 
a neighborhood of Q. 

Exactly the same analysis, with exactly the same result, can be applied to the 
second integral (over SR) in ( 13.131) . We have 

q>(uo,vo\u,v) ~ — A"log|a + v| as — e = u + v — >0— 

for some K. Thus in this limit there are three terms in ( 13.13b that will pick up 
a log(« + v) factor, the first term, and a contribution from the lower end of both 
integrals. The magnitude of these terms is determined by the behaviour of the 
characteristic initial data and, for generic data, cancellation will not occur. 

Based on the above observations we can write the principal part of a general 
solution 

Xj/(u,v) =T(u,v) + T_(u,v) + T(u,v) 
T(u,v) =(p(u ,vo;u,v) ^(«o,vo), 

ru 

T{u,v)= (p(u',v ;u,v)B[\i/](u',v )du', (5.8) 
T(u,v) = / (p(uq,v';u,v)B[y]{uo,v'] dv', 

Jvq 



with 



and 



B[y](«',vo) :=y u (u' ,v ) + h(u' + v ) V("',vo) 



B[v](woV] :=\i/ v (u ,v') + l(u +v') V(«o,v'). 



(5.9) 



, „ , , , 1/2 fu'+v\ l/2 ~f( u '+V)(u+vy 
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where for simplicity we have set F(y) := F(|, |; 1, 1 — y). 
We find 

t{ U ,v) = ¥ (u ,vo) ( u -^) 1/2 ( u -^) 1/2 f ( ; m ° +v ;; (m+v ; 

J \uq + vJ \u + vqJ \(uq+v)(vq + u) 
1 ("o+vq) v("o,vo 



n ( M0 _ M )i/2( M + V0 )i/2 
-:T{u) log \u + v\ 



log|a + v| 



and, for the first integral term, 

T(u,v) 

J»o ; V"-" + £ / \u + vq J \(u ' -u + e)(v Q + u) J 

7T y„ \U'-U J \ U + V / 

=: log|« + v|. 

The calculation for r(w, v) is similar. 

This establishes that the leading term in the asymptotic expansion of the func- 
tion \jf is of the form 

y(w,v) ~ f(w) log|M + v| as-e = a + v-^0, (5.10) 

and similarly for the solution c 

c(u,v) ~C(w) log|a + v| as-e = a + v— ^ 0, (5.11) 

where C and W are certain functions that are given explicitly by the above formulas 
and are generically non-vanishing. These coefficients make sense as distributions: 
for instance T_(u) is defined the limit of the sequence 

7»:= (\w]iu', V0 )(4^\ l2 ( U ^\ l2 d u'. (5.12) 
Ju \u' -u + e J \ U + Vq ) 

Each term Tf ( • ) belong to L 2 (as follows from the regularity assumed on the data) 
and, as a sequence, converges to the distribution T_(-). Hence, for both C and W 
we have the existence of C e , W 8 such that 

(C,W) = lim(C e ,Y e )(-), (C £ ,r)(.)eL 2 , (5.13) 

£->0 



where the convergence holds in the distributional sense. 
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5.4 Proof of the main result 
Blow-up analysis 



We may infer the leading term in the asymptotic expansion of the metric function 
a from ( 13.3b as 



a(u,v) ~A(«) log|M + v|, A = C + 
More precisely, one has 

A = limA £ , A E e L l 

£->0 



(5.14) 



(5.15) 



where the convergence holds in the distributional sense. The behavior of the metric 
function b is given by ( 13.21) as 



ilog|«- 



(5.16) 



and so we can construct the leading terms in the asymptotic expansion of the 
metric 

g ~ |« + v^dudv - |w + v\ l+2C dy 2 + \u + v\ l - 2C dz 2 ) . (5.17) 

The Ricci and Weyl curvature tensors are best described using the Newman- 
Penrose formalism and explicit formulae for the NP scalars <P mn , A and W n were 
given in Section |4] After some calculations, the Ricci and Weyl invariants as \u + 



are given by 



/? = 24A~i(4C 2 -4A-l)| M - 



-2-2A 



\R a pR al3 1 1/2 = 2a/2 I #00*22 + 18A 2 + 2<P n 2 - 4| 4>, 



1 

4\/2 



4A - 4C Z 



'oil - 

i|(i + | M +vr 4A ) 1/2 | M - 



1*02 
,.1-2 



2 1 1/2 



and 

\r n rr a P"f 5 \ l l 2 

\ K -apyo'^ 



= 4 1 3W + %V 4 - 4«P, W 3 
r(2A+4C 2 - l) 2 |a + 



1/2 



v\- 4A + 4A 2 C 2 \u- 



4.4 



1/2 



|u+v| 



-2 



We emphasize that #00, #01 , #11 , #22 belong to L 2 so that the scalar (R a pR a ^) ^ 2 
belong to L 1 . On the other hand, the Ricci scalar has an expansion whose the 
principal terms's coefficient is the limit of functions in L . 

A necessary condition to avoid the blowing up of the curvature of the space- 
time at u + v = is that these leading terms vanish, which occurs if and only if 



C / -i=A = 0. 



(5.18) 



This condition implies W = 0. Even if these very special conditions are fulfilled 
there is no guarantee that spacetime would not blow-up — one would have to in- 
vestigate the next order terms. We conclude that in the case when the above coef- 
ficients do not vanish, u + v = is a strong curvature singularity. 
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Finally, we emphasize that by perturbation of the data in the natural energy 
norm (corresponding to their assumed regularity), we can always ensure that the 
above coefficient do not vanish. Indeed, if ( 15.181) holds, then by adding a constant 
a to the characteristic initial data for iff, we see that 

B[Y+a](u',v ) :=B[^](M / ,vo) + i(a / + v )" 1 a, 

so that the coefficient on the singularity (given by the limit of ( 15.121 )) is changing 
accordingly by a (non- vanishing) term proportional to a. In turn, C 2 — j and A 
must be non-vanishing for all sufficiently small a, at least. This completes the 
proof of Theorem [2] 



5.5 Special solutions 

It is of interest to consider the special case of trivial initial data, that is, when 
there are no gravitational waves and no matter content. For simplicity let us set 
«0 = vo = 0. As noted earlier, r = e b satisfies dud^r 2 ) = nd, after normalization, 

r 2 = f{u) 2 +g{vf - 1 with f(u) = r(«,0), g(v) = r(0,v), and (/(0) = g(0) = 
r(0,0) = 1. 

When the free data vanish, we can solve the relevant Einstein equations and 
obtain that / and g vary linearly, that is, /(h) = 1 + (Xu and g(v) = 1 + j3v for 
some constants CC,f}. We distinguish between two cases whether these constants 
have the same sign or opposite signs and we make here specific sign choices since 
the other cases can be deduced by time or space reversals. In addition, by taking 
advantage of the transformation (w,v) h-> (pu,v/p) we can finally consider Case 
1 : a = P > and Case 2: -a = J3 > 0. 

With vanishing free data, the metric reads g = e 2a (du 2 + dv 2 ) — r 2 (dy 2 +dz 2 ), 
while the contraint ( 14.101) takes the form f^> = e~ 2a r~ 2 d u rd v r and the Bianchi 
identities become <9 u (r 3 f^) = <9 v ,(r 3 f<>) = 0. Hence, for some constant m we have 
r 3l F 2 = -2m. 

Consider first Case 1 above. Along ,jY_ one has ¥^ = r~ 2 d u rd v r and 

r=l+au, rd u r = a(l + au) , rd v r = a. 

Therefore, r 3 p = — a 2 and 2m = a 2 > and we obtain e 2a = (1 + ecu) (I + Ctv) /r. 
In this case, we have u > —1/a and v > — I /a. Coordinate singularities arise 
at u = — I /a where d u r vanishes, as well as at v = —1/a where d v r vanishes. 
More importantly, a genuine (geometric) singularity arises at r = 0, that is, on the 
"quarter circle" determined by the conditions 

{u + \/af + {v + \/a) 2 = l/a 2 , u>-l/a, v>~l/a. 

If now we make a conformal transformation of the metric we can set a = 1 
and apply the transformation u i— > u — 1 and v i— > v — 1 , so that the two metric 
components are given simply by 



e 2a = -, r 2 = u 2 + v 2 -\. 
r 
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The quotient manifold Q then corresponds to the part of the positive quadrant 
which lies outside the unit circle. The maximal extension is obtained by introduc- 
ing U = u 2 /2 and V = v 2 /2 which removes the coordinate singularity. By setting 
U = T — x and V = T + x the metric becomes 

g = (4t- iy l/2 (dT 2 -dx 2 ) - {Ax-\){dy 2 +dz 2 ), 

and the extension covers the domain T > 1/4. Finally, setting 3f = (4t — l) 3 / 4 and 
making the transformation (x,y,z) (3 l / 3 x, 3~ 2 / 3 y, 3~ 2,/3 z), we see tha the metric 
takes the well-known Kasner form 

g = dt 2 - r 2 l 3 dx 2 - t A l\dy 2 + dz 2 ). 

In addition to the assumed symmetry made for general solutions, this special case 
enjoys a rotational symmetry in the (y,z)-plane and, furthermore, the Kasner met- 
ric is easily checked to be conformally invariant so that the constant a has been 
absorbed and we arrive to only one geometrically distinct solution rather than a 
one-parameter family of solution as one could have expected. 

Concerning Case 2 mentioned earlier, a very similar analysis can be performed 
which leads us to the metric expression 

g = x- 2 'W - r^dx 2 - x 4 / 3 (dy 2 + dz 2 ). 

This is now a static solution with a timelike singular boundary. Again, this metric 
is invariant under homotheties, so that the constant j3 is absorbed and there is again 
only one geometrically distinct solution. 



5.6 An alternative approach 

Another approach to studying the behavior of solutions on the coordinate singular- 
ity is provided by imposing data directly on the singularity hypersurface u + v = 0, 
as we now discuss. Following Hauser and Ernst 1201 . we start from the general 
form of an EPD equation 

(a + v)w + ay/„ + j3yA, = 0, (5.19) 

where a, j3 € (0, 1) are constants. It is easy to check that (it — a)~^(v + o)~ a is a 
solution for any constant a. Then, superposing such solutions we see that 

Jx (w- ctj^v + ct)™ 

is also a solution, where A = A (a) is an arbitrary function and X and Y are con- 
stants. A formula due to Poisson allows also to take X = u and Y = —v. Hauser 
and Ernst looked at the case X = uq, Y = u with a = j3 = I, and observed that 



A(a)da f A(a)^/vo + o do 

is a solution of ( 15.19b with a = j5 = A. 



. . /•» A{a)da f" v 

y/(u,v) = / ^— = / r — ~ (5.21) 

Juq V" — 0V V + tT Juq \/U — Oy/V + a 
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For the moment we accept this assertion, and evaluate the solution when v ■■ 

Vn, 



f u A(a)da 

y(u,v ) = / V— ■ (5.22) 
Juq V u — ° 

Now the range of this integral is the line SQ, and the left hand side is (that half of) 
the characteristic initial data which we wish to impose on SQ. Note that we are 
assuming \ff(uo, vn) = 0. If iff represents the velocity potential this is not a prob- 
lem. If \ff represents the metric coefficient c we may need to scale the ignorable 
coordinates y and z to achieve this. 

As Hauser and Ernst pointed out, ( 15.22b is an Abel integral equation which can 
be solved explicitly for A (a): 

A{o)= l *n^ U . (5.23) 

Since a = j5, we can repeat the argument by interchanging u and v and this leads 
us to the formula 

. , r v B(a)JuQ + a da - , Id r° \ir(u ,v)dv 
y (u,v = / — J^Ll , 5(a) = -— / yv A— , (5.24) 

■Jv \Jv — 0\Ju + o naaJvQ \Jo-v 

produces a solution of the EPD equation which fits the characteristic initial data 
on SR. By linearity, the general solution to the characteristic initial value problem 
is then given by summing (15.21b and (15.24b . We thus have obtained explicitly the 
dependence of the solution on the initial data. 

We do need to check that ( 15.211) and ( 15.24b are in fact solutions which can 
be continued up to the singular line u + v = 0. (This is not a problem for the 
Poisson-Appell version (15.20b .) Concentrating on (15.21b we change the variable 
of integration so that the integral is over a fixed range by setting 

a = (1 — t)mo + tm, 

f = U-U , /l = l-T, f 2 =V + Uo + T(u-Uo). 

Then Hauser and Ernst assert that 

y(u, v) =r^^ (5.25) 

JO V/lV/2 



is a solution for arbitrary A (a). Suppose we define 

G ^ = A -WrrFTr- (5 - 26) 

V/0ilV/U2V/2 

Then it is straightforward to verify that ( 15.25b is indeed a solution for arbitrary 
A (a) provided the "surface term" 

-1 T=l 

G(t,m,v) (5.27) 

J T=0 



The characteristic problem for plane symmetric spacetimes with weak regularity 



35 



vanishes. This is a delicate matter. Near T = 



A{uq)^/u-uqX 
G(t,w,v) — , 

and this vanishes as t 4- 0+ provided v + uq^O. Near t = 1 



G T,M,V) - - 



(a + v)V« + v 

Away from the singular line m + v < £ < 0, the limit as T f 1 — is straightforward, 
and ( 15.21b does indeed furnish a solution of the EPD equation. However if we want 
to evaluate \ff(u,v) near the singular line, (u + v) f 0—, the (somewhat ill-defined) 
order in which we take the limits is crucial. 

Suppose we accept that ( 15.211) is a solution. In general we cannot calculate 
its value. As a special case suppose that A(ct) is a constant, say ia. (The i factor 
allows for the fact that v + a < in our scenario.) The indefinite integral is easily 

seen to be 

alog(V"— ffy— v— a + 2<7 — u + v), 

and 

!//■(«, v) =a(log(a + v) — log(2\/ a — «oV — v — M o + 2wo — m + v)) , (5.28) 

which exhibits the log(a + v) singularity that we expect. We conclude with the 
following question: The condition for avoiding the log singularity is that A (a) — > 
sufficiently fast as a f u—. Clearly A (a) must behave in the same way. It would 
be interesting to use ( 15.23b and relate this to the behavior of \jf(u,Vo). 

Let us make a final remark when analyticity of the data on the singularity is 
assumed. We recast the EPD equation ( 13.5b in the form 

L[y]=M[y], (5.29) 

where L[\ff] = \\r tt + \j/,/t and M[\jf} = i#a. e.g., L[t n ] = n 2 t n ~ 2 , M[x"} = n(n - 
l)x n ~ 2 . Suppose we have a sequence of functions {...,(«„(?),...} with the prop- 
erty L[(O n ] = (0„-i, and a sequence of functions {. . . ,(j> n (x), . . .} with the property 
M [</>„] = (j)„+\. Now, ignoring all questions of convergence, it is easy to see that 
y/(t, x) — Em=-°° <u "i(f) < /'"i( x ) satisfies the EPD equation. 

On physical grounds we do not expect iff to contain arbitrarily large negative 
powers of / as t — > and so we require (0„(t) =0 for all n < hq. Since m in the 
sum above is only defined up to an arbitrary additive constant we may set tiQ = 0, 
i.e., 

Xj/(t,x) = £ ft)„ (f )</>«(*) ■ 
«=o 

Clearly L[wo] = which implies (Oo(t) = A +5 log?, where A and B are constants. 
This in turn demands 

(O n {t) =A(a„t 2 ")+B(a n t 2n logt +b„t 2n ), 

where a n = 2~ 2n (n\)~ 2 and b„ is defined by a recurrence relation 

b n = (2n)~ 2 {b„-\ -4na n ) 
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with bo = 0. Suppose that at t = the solution of the EPD equation is 



Xff(t,x) =Aa{x) +S/3(x)logf. 



(5.30) 



Then, within the class of real analytic solutions one finds 



\jf{t,x) =A£a„f 2 ' , a„(x)+B£(fl„logf+Z?„)f 2 "j3„(x), 



(5.31) 



o o 



where 



Unix) = hr a W' AiM = hr 





For instance, the exact solution cited in Remark [T] is generated by A = 8, B = 0, 



Clearly 1// is regular at t = if and only if B = 0, and such a solution will 
give a regular iff and i// r on another Cauchy surface, say t = to 7^ 0. However any 
slight change to the data at t = to means that B changes away from zero, and so the 
resulting solution will become singular as / — > 0. We deduce that generic solutions 
of a Cauchy problem blow up as t — > 0. 



6 Propagation of curvature singularities 

6. 1 Choice of tetrad and general expression of the metric 

In this section we determine which kind of curvature singularity is allowed on a 
null hypersurface, denoted below by JVo, for solutions of the Einstein equations. 
No symmetry assumption is made in our general discussion. However, in each 
statement below, we specify what simplification is achieved for plane symmetric 
spacetimes. Hence, this section provides us a geometric derivation of jump rela- 
tions satisfied by the spacetimes constructed in the rest of this paper, which could 
also be established from the distributional formulation of the Einstein equations 
that we have discussed. 

A main source of difficulty comes from the fact that the metric induced on 
a null hypersurface (from the spacetime metric) is degenerate, with signature 
(0,-1,-1). In consequence, several fundamental concepts (linear connection, 
etc) must be re-visited. This was done first by Stellmacher [33], but we adopt here 
the more geometrical approach due to Penrose [31]. He pointed out that it is natu- 
ral to distinguish between three distinct types of geometries on a null hypersurface 
and, in consequence, to impose the condition that these geometries coincide when 
taking the limit of the spacetime metric from either side of the hypersurface. 

We may distinguish between the following three classes of spacetimes j% = 
U singular across a null hypersurface jYq : = n : 

• I-Geometry : the induced, degenerate metric only is continuous across jVo- 
This is the most general situation where most curvature scalars contain Dirac 
measure terms. 

• II-Geometry : a concept of parallel transport, along each integral curve 7, of 
tangent vectors to 7 is assumed to coincide on both sides of the hypersurface. 



and ot(x) =x 2 . 
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• Ill-Geometry : a concept of parallel transport, along each integral curve y, of 
tangent vectors to the hypersurface jV is now assumed. 

The original presentation by Penrose was based on the use of 2-component 
spinors. Our presentation will use the more familiar Newman-Penrose formalism 
based on a null tetrad, which avoids the explicit use of spinors, but restricts our 
consideration to II- and Ill-geometries. We will construct the spacetime by start- 
ing from an given arbitrary spacelike two-surface and considering hypersurfaces 
generated by the two families of null geodesies orthogonal to that two-surface. For 
each singular geometry, we will then investigate which singularities are admissible 
along these hypersurfaces. 

For the moment we assume that all of the quantities under consideration are 
sufficiently regular. Let be an arbitrary spacelike two-surface with intrinsic 
coordinates (x A ) with A = 1,2. Through each point P ,o £ ^0.0 there pass two null 
geodesies orthogonal to the tangent space T^q q. As P ,o varies (and provided 
they do not recross), these geodesies trace out two null hypersurfaces which we 
denote by ,jYq and JV^, respectively. Let u be a parameter along the generators of 
jV§, normalized so that 



The parameter v is similarly denned on jVq, and neither parameter need be affine. 

Let £P u fi be the (necessarily spacelike) two-surface determined by the con- 
dition u =constant in JVq. Through each point P u q € £? u a there passes a second 
null geodesic orthogonal to T g? u $ and (assuming again that no conjugate point 
arises), these geodesies trace out another null hypersurface JY„. We extend the 
definition of the parameter u by requiring that u remains constant on JYJ. In a 
similar way we construct a family of null hypersurfaces ,JV V and, now, u and v are 
two of the four coordinates on the spacetime Jt which, without loss of generality, 
we assume to be entirely covered by our construction. The two-surfaces on which 
u and v remain constants is denoted by ^ u ,v 

We extend the definition of (ji" 4 ) from ^0,0 to 3^ u a by requiring that ^ be 
constant along the generators of ,yfo. We then extend the definition of (Y 1 ) to the 
surface ^ HiV by requiring that ^ is constant along the generators of ,_A^- Note 
that, in general, ^ is not constant along the generators of Jf v , the exception being 
the initial hypersurface 

Importantly, we can carry out these constructions both for v > and for v < 0, 
and these are (locally) distinct and lead to two pieces of our spacetime, denoted by 
and M~ , respectively. By construction, the coordinate charts are continuous 
at,Ao = n but no further regularity is available in general. 

Consider, for instance, the manifold (with boundary) . The surfaces JV^ n 
Jt + are null hypersurfaces in on which we have the null covector field 



and the vector field n a is parallel to the generators of ,yVj , along which only v 
varies. Thus, there exists a scalar function Q > such that 



u = on &> K 



'0,0- 



n a dx a = du, or equivalently n a 



= V a u in^# + , 



A := 




— in Jl . 



dv 



(6.1) 
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Similarly, the surfaces jY v n Jl + are null hypersurfaces on which we can de- 
fine the null covector field 

l a dx a = Q~ l dv in^ + , 

which enforces the condition l a n a = 1 in ^# + . The vector l a is tangent to the 
generators of ,jV v along which u and x A change. Enforcing n a l a = 1 we deduce 
that 

D:=ICC ^ = ^-+ cA 4a in ^ + ' ( 6 - 2 ) 
ox a ou ax A 

where C A are two real-valued, scalar functions. 

Finally, consider the tetrad vectors m a and m a . Since n a m a = l a m a = by 
definition, they cannot have u and v-components, and we thus set 

8:=m a ^=P A ^ in^+, (6.3) 

where the P A are two complex scalars. Note that we are always allowed a spin 
transformation of the form 

m a ^m a = e ie m a , m a h+ M = e~ ie m a , 

where 6 is any real-valued scalar field. This corresponds to the transformation 
P A n> P a = e' e P A , and thus P A and P^ contain three real degrees of freedom, 
only. So, we define P A and P^ by the conditions 

P A P A = PaP* = 0, PaP* = PaP* = 1. 

Then, defining 

m a dx a = (P B C B ) du - P A dx A in Jt + 

completes the tetrad since m a l a = m a n a — m a m a — and m a m a = — 1. 

Finally, we can obtain the general expression for the spacetime metric on ^# + , 
gab = 2Z( a n fo ) - 2m( a 7%j, and summarize our conclusions, as follows. 

Lemma 1 In the coordinates (u,v,x l ,x 2 ) constructed above, the metric g in both 
parts takes the general form 

g=-2\P A C A \ 2 du 2 + Q~ 1 dudv - h AB C A dudx 8 - h AB dx A dx B , 

where 

hAB = 2P(aPb) 

is a real-valued metric on the two-surfaces £? u ^ v . 

When plane symmetry is imposed, one can choose C A = identically, and the 
metric only depends upon the function Q, and the symmetric 2-tensor hAB 

g = Q~ x dudv — hAB dx A dx B ; 

moreover, all 8- and 8 -derivatives vanish. 
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Note that the metric has six explicit degrees of freedom. In this sense the co- 
ordinate choice is optimal. Also on ^% we have C A = by construction. 

Penrose's classification of the geometries allowed on jY§ can now be described 
as follows: 

• I-Geometry : The function Q to be discontinuous across jVq, while the coeffi- 
cients C A = and P A are continuous. 

• II-Geometry : The function Q is continuous. 

• Ill-Geometry : The parameter v is affine, so that Q is constant on jYq and 
without loss of generality one can assume that Q = 1 on 



6.2 Jump relations for the NP connection scalars 

We now discuss solutions that may be singular across the hypersurface jYq, which 
has v = by definition. We will need to handle first- and second-order derivatives 
(in the variable v) of the coefficients C A and P A and, for instance, our calculation 
will involve the jump [C A v ] of the derivative of coefficient C A across the hyper- 
surface jV§. 

With Penrose's standard notation [28 34 1 and using the expression of the met- 
ric in LemmaQ] we obtain 

a = {P^ - \P A 8P^ - \8\ogQ- \P A AC A , 
p = ^SP 4 - \P A 8T" - \8\ogQ-\P A AC A , 
y=\P A AP A -\P A AF A , 

e = \P A DP A + {PaDP^ - \P A 8C A + \P A 8C A - \D\ogQ, 

)c = 0, k = -P A AP A , 

ii = - 1 1 P A AT'-rp A AP A , v = 0, 
K = ~ l {8\ogQ-\P A AC A , 

and 

p = iPA/>^ + jPaDP* - \P A 8C A - ±P A SC A , 
a = P A DP A - P A 8C A , 
T=\8\ogQ-\P A AC A , 

where the derivatives D, A , <5, 8 are understood in the sense of distributions. 
We have therefore established the following result. 

Proposition 3 (Jump relations for Newman-Penrose connection scalars) With 
the notation above, the NP connection scalars K and V vanish identically while 
p, o, and e are continuous across the hypersurface The other NP connection 
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scalars may be discontinuous across with jumps given by 
[a] = -\P A [C% [j8] = [a], 

[y] = \PA[P?]-\PA[tl [X]=-P A \P$], 

[T]=2[S]. 

Wzen plane symmetry is imposed one has also a = f} = n = T = 0, so that the 
remaining jump conditions are 

[X] = -P A [t], 

[li} = -\P A [P?]-\P A ^\. 



6.3 Jump relations for the NP curvature scalars 

One more piece of notation is now needed. If a scalar / takes values on jV§ 
when approached from j& , we set 

while we have already defined [/] = / + — /~. We note the classical identity 

\fg]=f[g} + [f}g- 

o 

Recall that if / contains a Dirac measure term, its coefficient is denoted by /. 

We now impose Einstein's field equations in the NP notation, as stated for 
instance in P4l (p. 217). By the equation (a) (in the notation therein) we find 

<J> 00 = Dp - p(p + ee) - oo : 

whose right-hand side (according to Proposition [3]) is continuous across jVq, so 
that [<f>oo] = 0. Similarly, by equation (b) in [ 34 1 we find 

% = Do (2p + 3e-e)(7 

which implies [Wo] = 0. The coefficient <E>22 is determined by equation (n): 

<f>22 = -Aji - ji 2 - XX, 

o 

which with Proposition [3] implies that <J>22= — [m]- Similarly, equation (j) yields 

W 4 = -AX-2(n + y)X, 

thus Wa= -[X]. 
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The coefficient 4>o2 occurs in equations (g) and (p). Equation (p) involves the 
term Aa = <7 y , which might have a non-zero jump across Jfo (which could also 

o 

be computed explicitly). However, by equation (g) we find <£>2o= and 

[<f> 20 ] = D[X]- 28[a] - (p - 3e + e) [X] - a\p] - 7.(1% - a + j8 ) [a] . 

Next, W\ and 4> i are involved in equations (c), (d), (e), and (k). Equation (e) 
implies 

Wi=0, 

[fi] = D[a] - 3a[a] - (p - e + e)[5], 

while (d) gives 

<lio=0, 

[<P W ] =D[a]-(3p-e + e)[a]-o[a\. 

Equations (c) and (d) then provide no new information. 

The terms W 2 , <£ii, and A occur in (f), (h), (1), and (q). Unfortunately, (f) 
involves 8e while (q) involves Ap, both of which are tedious to compute. We can 
however deduce 

o o o 

W 2 =<Pn=A=0, 

and the jump relations 

pft] + 2[A] =D[n] - 28 [a] - (p - e - e) [p] 

-a[X]-2(%-^+ j8) [a] - 25? [a] , 

and 

-pft] + [A} + [<P n ] =8[a] - 8[a] -p[p] + cr[A] 

-(B+3j3)[a]-(3a + ji?)[5]. 
Finally, and <J>i 2 occur in (i), (m), (o), and (r). Equation (r) provides us 

o o 

with W 3= — [ot], while (o) gives <Pu= - [a] . Equation (i) is consistent with these 
results and gives no further information, but (m) reveals that 

-[%] + [*2i] = S[X]-S\pL] -(a+^M- (5 -3/3) [A] -2p[a] +2X[a]. 

We summarize our conclusions as follows. 

Theorem 3 (Jump relations for Newman-Penrose curvature scalars) Assum- 
ing that the hypersurface jY§ has type II-geometry, then the Dirac measure terms 
on it satisfy the following conditions: 

*P 4 =-[X], «P 3 =-[a], J>22=-[p], *i2=-[a], 

o o o o o 

<pQo=<Pw=<P2o=<Pn=A=0 1 

o o o 

«f / 2=f / l=f / 0=0. 
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If one strengthens the assumption to that of type Ill-geometry, i.e., [a] = 0, then 

o o 

in addition W3=<Pu= 0. When plane-symmetry is imposed the latter condition 
are automatically satisfied and the geometry is always of type III. If one consid- 
ers "pure gravitational radiation " which, following Penrose, requires no Dirac 
measure terms in the trace-free Ricci tensor, then one has [jii] =0 as well. 

These results can be interpreted in terms of physical quantities, as follows. For 
instance, a concentration (Dirac mass structure) in the energy density of the fluid 
requires a jump in the Ricci curvature, their relative strength being proportional. 
This statement, in particular, is valid as one approaches the future boundary of the 
maximal development (considered in the previous section) and provides a corre- 
sponding relation between the curvature and matter content of the spacetime. 
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